Fluctuation conductivity in superconductors 
in strong electric fields 



Todor Mishonov,"^^* Anna Posazhennikova" and Joseph Indekeu" 
Laboratorium voor Vaste-Stoffysica en Magnetisme, Katholieke Universttett Leuven, 
Celestijnenlaan 200 D, B-3001 Leuven, Belgium 
''Department of Theoretical Physics, Faculty of Physics, Sofia University St. Kliment Ohridski, 
, 5 J. Bourchier Blvd., 1164 Sofia, Bulgaria 

We study the effect of a strong electric field on the fluctuation conductivity within the time- 
f~| ' dependent Ginzburg-Landau theory for the case of arbitrary dimension. Our results are based on 

I the analytical derivation of the velocity distribution law for the fluctuation Cooper pairs, from the 

Boltzmann equation. Special attention is drawn to the case of small nonlinearity of conductivity, 
QQ ' which can be investigated experimentally. We obtain a general relation between the nonlinear con- 

ductivity and the temperature derivative of the linear Aslamazov-Larkin conductivity, applicable 
to any superconductor. For the important case of layered superconductors we derive an analogous 
relation between the small nonlinear correction for the conductivity and the fluctuational magneto- 
conductivity. On the basis of these relations we provide new experimental methods for determining 
both the lifetime constant of metastable Cooper pairs above Tc and the coherence length. A sys- 
tematic investigation of the 3rd harmonic of the electric field generated by a harmonic current can 
serve as an alternative method for the examination of the metastable Cooper-pair relaxation time. 
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Fluctuation phenomena in superconductors have been intensively studied for more than 3 decades since the discovery 
of the fluctuation smearing of the superconducting transition. The interest in fluctuations was resurrected when 
the high-tcmpcrature superconductors (HTSC) were found. In conventional BCS-like superconductors the transition 
temperature Tc marks a sharp dividing point between distinct regions of 'superconducting' and 'normal' behavior, and 
critical fluctuations are almost unobservable. Largely due to their high temperature, small coherence length and quasi- 
two-dimensional nature, HTSC show a signiflcant smearing of the critical transition and thus the effect of fluctuations 
^-H in the critical region of copper-oxides is much more pronounced and experimentally accessible. The conductivity, the 
specific heat, the diamagnetic susceptibility, the thermopower etc. have been observed to increase considerably in the 
normal state in the vicinity of the transition temperature in cuprate compounds. Importantly, for high-quality HTSC 
. samples the interactions between fluctuations are negligibly small, fluctuations can be described as Gaussian ones and 
■ their theory is rather simple. As a result nowadays the study of fluctuation phenomena in superconductors occupies 
an essential and significant part of the whole physics of superconductivity. For a contemporary review we refer to 
the work of Larkin aud Varlamovtl; cf. also the classical review by Skocpol and Tinkhama and the recent review by 
Mishonov and Penevcl especially devoted to Gaussian fluctuations in layered superconductors. 

In this paper we are concerned with the effect of fluctuations of the superconducting order parameter on both 
(~| the linear and the nonlinear conductivity near and above the bulk superconducting critical temperature. In order to 
O ■ account for the fluctuation conductivity (or paraconductivity) one needs to apply the time-dependent generalization 
^ of the Ginzburg-Landau theory (TDGL). The theoretical study of fluctuation conductivity dates from the paper of 
^ , Aslamazov and Larkin in 1968 based on the microscopic approach. Although a long time has elapsed, nonlinear effects 
have not been sufficiently investigated and the theoretical results obtained so far require specification. 

The aim of the present paper is to study the kinetics of the superconducting order parameter in a strong electric field 
• within a simple phenomenological approach and to provide new experimental methods for determining such important 
" ■ ' material constants of superconductors as the lifetime of fluctuation Cooper pairs tq and the coherence length ^(0). 

We adopt the ideology of TDGL theory and show that the velocity distribution law of fluctuation Cooper pairs can 
be derived from a Boltzmann equation. The cases of superconducting wires, thin fllms, bulk superconductors and 
striped and layered superconducting systems are analyzed. We are particularly interested in the effects of nonlinearity 
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of conductivity, which can be observed experimentally. We show that the strong critical behavior in the vicinity of Tc 
of the third harmonic of the electric field generated by a harmonic current, can be used as an experimental technique 



for the determination of the metastable Cooper 
This paper is organized as follows. In Section M 



jair lifetime constant. 

we discuss the kinetic Boltzmann equation for fluctuation Cooper 
pairs in strong electric fields and present the momentum distribution law for Cooper pairs in the normal phase and 
the expression for the fluctuation current. Section III is mainly technical and devoted to dimensionless notations we 
introduce to make our calculation less complicated. The derivation of the time evolution of the distribution function 
using canonical variables is given in the Appendix. Section IV explains the paraconductivity in layered materials, 
such as for example HTSC. We provide the analytical derivation of the conductivity correction in the normal phase in 
a Lawrence-Doniach superconductor. In Section ^ the Aslamazov-Larkin conductivity dependent on the electric field 
is described in detail for one-, two- and three-dimensional superconductors. We study strong electric field effects on 
conductivity as well as weak electric field influences below the critical temperature. The interesting case of a striped 
superconductor, which is becoming very popular nowadays, and that of a thick film can be found in Section VI. The 
conductivity of such a superconductor depending on the thickness of the film, is obtained. The expressions for the 
nonlinear conductivity correction can be used for experimental data processing for extracting the lifetime constant 
of Cooper pairs, as is shown in Section VII. A new experimental technique for probi ng tq, based on suggested 
measurements of the third harmonic of the electric field, can be found in Section VIII . Finally, a discussion and 
conclusions are given in the last section. 



II. BOLTZMANN EQUATION AND FORMULA FOR THE CURRENT 

As was already mentioned in Section |, the Ginzburg-Landau theory serves as an adequate tool for describing the 
fluctuation conductivity phenomena in superconductors. The important fundamental constant of the TDGL theory 
is To, which is proportional to the lifetime of metastable Cooper pairs in the normal state: 

r(6) = ^, (2.1) 

where e={T-Tc)/Tc. 

Within the weak-coupling BSC theory in the case of negligible depairing mechanisms the temperature-independent 
constant To satisfies the relation: 

At present the experimental value for tq, obtained for the layered cuprate superconductors, is in good agreement with 
the BCS theoryd. From an experimental point of view it is convenient to introduce the dimensionless ratio: 

(BCS) ' ^'^■^^ 
Tq 

which just characterizes the deviations of the experimental value from that which is theoretically derived from the 
microscopic theory for the case of small coupling and negligible depairing. Before establishing the relation between 
experimental results on fluctuation conductivity and the lifetime constant To the important question to be clarified is 
what the momentum or velocity distribution of the metastable Cooper pair actually is. 

The charge of the Cooper pair is |e*| = 2 |e| and its effective mass m* is unambiguously experimentally accessible 
through many principally different methods: electrostatic charge modulation of the kinetic induptanceO, surface Hall 
effectQ, magnetoplasma waves and cyclotron resonances, Doppler effect of Cooper pair plasmonsa and, which as most 
reliable, determination of the thermodynamic equilibrium electric potential related to the Bernoulli potentialEJ. The 
velocity of a pair along the a;-direction reads 

v.^^ (2.4) 

and along the y and z directions the velocities are correspondingly Vy = py/m'^, and Vz = Pz/m*. Note that we allow 
for the presence, oij rnass anisotropy. 

It was showntllEj within the general TDGL theory that, surprisingly, the momentum distribution law for the 
metastable Cooper pair is described by the classical Boltzmann equation introduced in physics in 1876 long before 
the electron was discovered. Nowadays this equation is discussed in many well-known textbooks on kinetics and solid 
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state physicsS. For a-detailed description of the relaxation time approximation to the Boltzmann equation see also a 
well-known textbookd. The Boltzmann equation in our context is given by 



| + e*E.|)n(p,t) 



n{p,t) - Up _ np{t) ^ rPT 

Tp Tp To 



(2.5) 



This kinetic sanation for fluctuation Cooper pairs has been applied to the fluctuation Hall effect in thin supcrcon- 
ducting-filmaiil and to the fluc tua tion paraconductivity within the framework of time-dependent Ginzburg-Landau 
theoryllj, cf. also Ref.EI. In Eq.(2.5) np{t) — n{p, t), and we introduce the dimensionless number nx — ksT/aQ, where 



2miea{0) 2m*42(0) 2m* eM 
is in fact proportional to the first coefficient in the Ginzburg-Landau expansion 

a(e) = aoe 

and 

Up = 

Sp- ^I 



(2.6) 
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is a standard equilibrium distribution, where the energy spectrum of the layered superconductor with layer spacing s 
is given by 
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2m*, 
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(2.9) 



and the , chemical potential is defined through —p,{T) — a{e). The momentum-dependent relaxation time Tp 
edBE3b 



obtainecOEj by TDGL theory reads as 



Ep + a(e) ■ 



(2.10) 



Let us assume that the electric field is constant and applied along the x-direction in the a6-plane: E = {E^, 0, 0) 
const. In this case from the Boltzmann equation (^.5[) follows the stationary momentum distribution law 



i{px\e±,E^,T) = 



n'/2m*ea{0) Jo 
X exp 



exp ■ 



H/to fl \To 



h'/2mieai0) 

1 / e*E,To 
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where 



£A-{Py,Pz) 



2m*, 



1 — cos 



(2.11) 
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and t he physical meaning of the dummy parameter t as time is revealed in the Appendix. For zero field E^ = from 
( ^.11 ) we naturally obtain 



ksT/ao 



iUO)Px/fir +£±/ao + e 



(2.13) 



Once we have the momentum distribution law for the fiuctuation Cooper pair, we can immediately proceed to the 
expression for the stationary fiuctuation current for arbitrary dimension D: 



e v,-=y e v,n(p)^^, V, 



d^p dsp 
dp 



(2.14) 
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where V = L^LyLz is the volume of a system. For the current along the cc-direction in a layered superconductor we 
have 



dt 

— exp ■ 



3 \T^IU^)J \ro 
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Trh/s 



dpz f°° dpy 



2TTh ./ _ 2TTh 



X exp • 
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Px \ t f £.a{0)Px\^ , e*E,^(,a{0) f t\'^^a{Q)Px \ dp,. 



ml To 



^/ To V '^0 
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(2.15) 



A similar forpijila for the current was presented by Gor'kovlli and modified for layered superconductors by Varlamov 
and Reggian: 



Hi. 



III. DIMENSIONLESS VARIABLES 



To make our further calculations a bit less cumbersome we introduce a number of dimcnsionless variables: — 
£.a{0)px/fi and ky = £,biO)Py/^ G (—00,00) are dimcnsionless momenta, 9 = spz/h G (— tt, tt) is the Josephson phase 
and u = I/to G (0, 00) is in fact the renormalized time. We also introduce 
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ao 



2ec(0) . spz 
sm — — 

s 2% 



kl + Loie), 



(3.1) 



where for the well-known Lawrence Doniach model ijj{9) = Sz{pz)/ao = ^ cosO) = rsiv? |. The factor r = 
(2^c(0)/s)^ originates from the parametrization of the effective mass in c-direction, m*, for an anisotropic GL model. 
It is convenient to deal with some dimcnsionless variable /, proportional to the electric field 



e*ExTo 



and with the parameter (7, proportional to the nonlinear electric field correction 

P 1 fe*ExUOy^ 
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(3.2) 



(3.3) 



With the above notations we can now rewrite our expression (2.11) for the momentum distribution of the metastable 
Cooper pair in the layered superconductor: 



exp 



{kl + e + w)u + fk^v^ - -^fu^ 



du. 



This function is a solution of the static Boltzmann equation which can be written in the form (fc = kx 

^dn{k) 



f- 



dk 



{k^ + t + w) n{k) + TiT- 



(3.4) 



(3.5) 



Let us mention some details concerning the derivation of expression (3.4). Direct solution of the Boltzmann equation 
gives 



i(fc) = exp 



fc3/3+ (e + w)fc 
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TiT / exp 



no exp 



k^/3+{e + w)k \ dk 
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(3.6) 



^From here it is easy to see that the boundary condition is n(fco) = no. Then writing k — k — uf, introducing 
u = [k — k)/f and using 
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(3.7) 



we obtain from Eq.(3.6) 



n(fc) = riT 



ik-ko)/f 



-uq exp 



exp 

fc3/3+ {e + w)k 
1 
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du 



exp 



fcg/3 + (e + w)kQ 
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(3.8) 



In the limit ko — > — cx) for / > or fco ^ oo for / < we have n(—oo) = n{oo) — — tiq and we arrive at Eq.( p.4| 
Finally the expression for the current Eq. (2.15) in our notations for the layered superconductor takes the form 



Jx = e ^— ^a(O) 
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de 



X exp 



{- 



0)] u] 



X 



f°° dkx 
I kx exp { -ukl + fu^kx ] , ^ = crj,j;£;:r 
J-00 27r4a(0) 



(3.9) 



It is informative to compare this result with previously proposed expressiongli3o. Our further analysis of the fluctu- 
ation conductivity is based on this general formula. 



IV. PARACONDUCTIVITY IN A LAYERED METAL 



As we already pointed out in Section | the fluctuations should be more important in high-temperature materials 
due in part to their high transition temperatures and extremely short coherence lengths (on the order of 10 A). 
Thus the study of Gaussian fluctuation effects on the conductivity in layered compounds is indispensable and of 
primary importance. In this section we examine in detail the electric field influence on paraconductivity in layered 
superconductors and derive analytically the correcti on to the conductivity in the normal phase. 

The integration over kx in the current expression (3.9) can be easily performed: 



j kx exp {-ukl + fu'^kx} dkx = ^V7™/exp Q/^u^^ 



Hence we have for the current in the layered superconductor the expression 



Jx 



e*kBT 



de 



dk,. 



27ri/2?i^ 27rs 271^^(0) Jo 



du\ 



X exp 



-(1 - cos 61) 



u fu^ 

12-' 



^xx^x 



and, consequently, the fluctuation conductivity reads 

2 e^kBTToUiO) 



de 



X exp {- [e + ky +UJ 



dkii 



du\fu 



After integrating over ky we obtain from ( |4.3| ) 



^ / d-uexp{- [e 



+ uj{e)]u- gu^} ^jx/Ex 



Thus the conductivity in the Lawrence-Doniach (LD) model can be expressed in the way 
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du exp 
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(4.1) 
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(4.5) 
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or, more conveniently, 



e'fcBrTo^,(0) 
nhh^biO) Jo 



du exp {— (e + r/2) u — gu^^ Iq {ru/2) , 



where 



1 



Jo{ix) 



(4.6) 



(4.7) 



is the Bessel function for imaginary argument. 



In the case of zero electric field = the integral over u in Eq.(4^) is easy to compute since 

d9 1 1 



exp|— e+^(l — cos6') u| 



27r(e+0-§cos0 ^e(e + r) ' 



and we obtain the Lawrence-Doniach formula for the fluctuation conductivity of a layered superconductor 



o-LD(e;0 



(4.8) 



(4.9) 



Hereafter we define, for simplicity, cr(e;r) = iT(e, 0;r). ■— . 

It is simpler to deal with the dimensionless function (cf. Varlamov and Reggianitil Eq. (7)) 

(e, 5; r) = = V<^) r i - + " - 5^^'} ^0 (™/2) d". 



Then Eq. (4.6) reads 



c^LD (e. 5; r) = CTlo (e; ^^d (e> 5; 



(4.10) 



(4.11) 



It is important to note that the product of two exponents in ( 4.10 ) can be expressed in terms of the third-derivative 
operator acting on the function of e 



exp (— exp (— eu) = exp | 9'^~3 ) "^xp (— eu) . 



(4.12) 



This key observation immediately leads us to the simple general relation between the nonlinear fluctuation conductivity 
and the linear one in evanescent fleld, cr(e), 



CT(e >0,g) = exp ( g-^ ) cr(e). 



(4.13) 



which is similar to the relation for the magnetoconductivity, derived in Ref.i. For our dimensionless function (4.1C) 
we get 



1 



Q3 



9e3 



In the limit g — s- the electric-fleld-dependent conductivity can be written in a simple form 



a(e, /) = a(e) + Aa/(e, f)^i^l + g—j a(e), 
and thus the nonlinear electric-fleld correction to the conductivity reads 

f2 Q3 



(4.14) 



(4.15) 



(4.16) 



Now we can derive the conductivity, dependent on the electric field in the Lawr ence -Doniach model of a layered 
superconductor. In principle we may start either from the right-hand-side of Eq.(4.8) (i.e., we can take the third 
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derivative after averaging over the Josephson phase), or we can perform the differentiation before integrating over 9. 
The latter is a useful method for the evaluation of complicated integrals necessary for experimental data processing 
and leads to 
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1 93 1 



[e+r(l_cos6')]'' 6 de^ ^e{e + r) 
Finally for the Lawrence-Doniach conductivity we have 

where 



2 (e + r)3 



and other notations are 



\e*E,\U^)T{e) _ I/I _ |i?,| 



\e*mr{e) 



£3/2 

i?,(0)e3/2, EM 
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12 



|eie(0)ro^ 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



, . . Ca,fc,c(0) r /2ee(e) 



(4.21) 



Our dimcnsionless <j function can now be rewritten in the scaled variables, as can be seen from ( 1.10| ), 

■^LDle^S;^) =<rLD(l,5£;?'e). (4.22) 

In conclusion of this section we obtain the formula for the electric-field-dependent correction to the conductivity in 
the Lawrence-Doniach model for a layered superconductor: 



5 ^31 



7rrs6(0) 



E(e + r)] 
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(4.23) 



V. ASLAMAZOV-LARKIN CONDUCTIVITY FOR D-DIMENSIONAL SUPERCONDUCTORS 

In this section we are concerned with the derivation of the Aslamazov-Larkin electric-field-dependent conductivity 
for bulk superconductors (Z? = 3), thin films {D — 2) and wires (D=l). The results can be generalized for the case 
of arbitrary dimension. 

We start with a one-dimensional superconductor. In order to derive the expression for the one-dimensional fluctua- 
tion conductivity we should set Py = = in ( |2.9D and exclude from Eq. (4.3) the integration over the perpendicular 
component of the momenta 



, , / hky hO 

(the electric field is as usual parallel to the x-direction) . We immediately get the final expression 

2 e^ksTToUO) 



flD 



.1/2 



■ / -v/uGxp {— eu — gu^l = — ^. 
Jo Ex 



(5.1) 



(5.2) 



For the derivation of a two-dimensional conductivity, induced by the fluctuations, we have to cancel averaging with 
respect of the Josephson phase 9 in Eq. (4.4) and we find 



^20 = 



e^ksTraUO) 



poo 

/ exp {— eu — 51*3} du 
Jo 



,-(2D) 
Jx 

Ex 



(5.3) 
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The case of a layered superconductor turns into the bulk (3D) one in the limit r — > oo, wh en t he distance between 
the layers tends to zero. Taking into account the asymptotics for the Bessel function in Eq. (4^), cf. RefO: 



loix > 1) 



exp(~^^)/„(^u)«-^, 



we obtain the bulk fluctuation conductivity expression 



e^kBTTQ^O) 



Jq V" 



du j. 



(3D) 



(5.4) 



(5.5) 



These results bring us naturally to the general expression for the fluctuation conductivity for arbitrary dimension 
D and T above Tr : 



0-0(6, .g) 



e^fcBTToeg(O) 
2^-27r^/2?i2^i)(0) 



/•oo 

/ exp {— eu — gu'^} 
Jo 



du 



(D-2)/2 ' 



(cf. Dorseyta (1991), EqJ47), for e > 0). 

One can express Eq.( ^.6| ) in a form similar to Eq. ( 4.11 ): 

CTD(e,/) = CTD(e)<rD(ffe), 

where 

. (A 4r (^) e2 r(e)gg(.) 

is the fluctuation conductivity in zero electric fleld, and 

^a(e) for D = l, 
.D, .) U^Me) for i? = 2, 
^ ^ ^ - ^ Ca(e)6(e)Cc(e) for - 3, 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



The function 



'^D(ffe) 



^-'"^e-^-^''""dv, v = tu 



(5.10) 



is a dimensionless function containing the electric-field d epend ence of the fluctuation conductivity. 
For small electric field /e ^ 1, in accordance with Eq ( 4.15| ), we have: 



^D(/e) « 1 - 



A-DX (Q-D 



-D\fl 
12' 



(5.11) 



A. Strong electric field expansion 



In this subsection we study the conductivity analytically in the limit of electric fields that are large compared with 
the reduced temperature distance to the critical point. To this end we introduce the scaled variable e/ — 12^^^ ^/P^^ 
and focus on the regime of small e/. The sign of e/ is arbitrary. Our analysis applies to T < Tc as well as T > Tc. 
Changing the integration variable u in Eq. ( |5.6| ) to f'^^^u allows us to extract the dominant field dependence at the 
critical point. Subsequently we use the series expansion of exp(— eu) and integrate it term by term for obtaining a 
series in e/. We arrive at the result 

Me,g)^ I e^^{-eu-gu^]^jj^^=AL^^^^Y.—lr^[ ^ )• (5-12) 







71 — 
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The convergence is very good for |e/| < 1. For instance, for e/ = 0.2 the errors in zeroth, first, second and third order 
are, respectively, 12%, 1%,0.07% and 40 ppm (parts per million). This series leads to the following expression for the 
conductivity 



with 



and 



aD(e,/) = aD(/)SD(e/) (5.13) 



oo j p ( 2n+i-D \ 

^B(e,)^i+E ,r r(4^) ^'-''^ 



is the fluctuation conductivity at T^, cf.llHli. 



3(D+2)/62(4D-10)/37ri3/2;j2_^D (Q) 

(D+2)/3 

(5.15) 



B. Weak electric fields below Tc 



In order to study the fluctuation conductivity for temperatures slightly below Tc and in small electric fields we 
must take into account that the limit of zero field is singular due to the occurrence of bulk superconductivity. The 
contribution to the conductivity that we calculate in this subsection must be well separated physically from that due 
to the onset of bulk order. The approximation scheme we develop here is a good one for relatively weak electric fields, 
which satisfy at T < Tc the condition 

_ I/I 2|e£;.|ga(0)ro 

= ^ = n\e\y^ < 1- (5-16) 

Especially below Tc we have to take into account that the temperature distance to the critical point e is renormalized 
as can be seen in self-consistent mean-field-like approximations which we will briefly consider later. The origin of 
this effect lies in the non-linear character of the TDGL equation, which cannot be neglected for high densities of 
fluctuation Cooper pairs. The selfconsistent approximation decouples the non-linearity, resulting in a linear problem 
with a modified e, to be denoted by Cr- To alleviate the notation we will postpone this substitution until the end of 
this subsection. 



For T < Tc we write e = — |e| and to the integral in Eq. (5.12) we apply the Gaussian saddle point approximation 
for weak electric fields. This amounts to looking for the maximum of the argument of the exponential function, since 
the remaining factor is an algebraic function of u and therefore slowly varying. Defining 

F{u)^-\e\u+lj^ (5.17) 
the saddle-point approximation can be written as 



e--^")-e-K)^_|_,(._.„), (5.18) 



where uq is the minimum of F(u), given by 



uo = 2^; (5.19) 



the time interval = tqUq has a transparent physical meaning, |e*£'a;| tg — 2?i/^a(e). In order for this approximation 
to be accurate, the condition 
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F"{uo)ul > 1 (5.20) 

must be satisfied, which is equivalent to 

/. « 1. (5.21) 

This condition can be seen to arise from two requirements. Firstly, already in Z? = 2, in the absence of the power of 
u in the integral, the validity of the saddle-point approximation requires 

exp (^^^F"{uo)ui^ « 1 (5.22) 

in order for the integration interval to be extendable to (— cxd, oo). Secondly, in the presence of slowly varying additional 
factors in the integrand, we must check the consistency of the approximation by performing a Taylor expansion of 
the algebraic function of u about uq. If we denote this function by G(u), we may approximate this by the constant 
G(uo) in the integral, provided G{u) deviates only weakly from linearity in a neighbourhood of width wf around uq, 
where wp is the standard deviation of the Gaussian function. This is fulfilled when 

G"{uo)wj,/G{uo) < 1. (5.23) 

Since wj^ = l/F"{uo) and G{u) is simply a power of u, this condition coincides with Eq. ( 5.2C| ). 

Within the range of validity of the saddle-point approximation we thus arrive at the following result for the integral 

/ o \ (3~D)/2 / 4 X 

JD{-\e\J)^V^\e\-^'-''^/'^jJ ^^P(3^)' e2//^»l. (5.24) 



According to Eq. (5.6) and taking into account the renormalization of e to we obtain 



The replacement e ^ derives from the use of a Maxwell-type selfconsistent approach (cf. Ref.i) for solving the 
Boltzmann equation. The implicit equation relating the renormalized and bare parameters is given by 

Mo r *\ ^nM2 f d^P 

with 



* Le*A„fc(0)]^ / -_^n(p,E,e,), (5.26) 



d 



(5.27) 



where details of the procedure of ultraviolet regularization will be presented elsewhere. 

The most significant aspect of our result is the dramatic exponential increase in the fluctuation conductivity for 
small electric fields. Therefore, for <C 1, the fluctuation part becomes of the same order of magnitude as the 
normal-state background, fj£){—\e\,f) ~ ctat, and the fluctuation conductivity will no longer be just a perturbation 
but a significant part of the total conductivity 

CTtot = (TD{e, /) + <7Ar(r) = jtot/E,. (5.28) 

This extraordinary increase of the conductivity naiurally leads to a minimum in the current as a function of the 
applied field, in agreement with the Gor'kov analysisEj that the current-voltage characteristic must have a perfectly 
noticeable section corresponding to negative differential conductivity. At a suitable potential difference between the 
ends of the film, the generation of radiation should be observed. Indeed, for every e < the generation of radiation 
will start for electric fields lower than the critical one, i?gcn(£ = ~kl)i determined by the criterion 



djtotiEx) 



dE^ 

which is equivalent to 



= 0, (5.29) 
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/gen ^'^(-I'^l'/) 



+ a(-|e|,/gen) = -aAr(T). (5.30) 

/gen 



Within the selfconsistent approximation Eq. (5.27) it is easy to obtain theoretical formulae for the case when the 
fluctuations are nonlinear 

EgonH^rl) = i?c(0)/gc„(-|e.|), (5.31) 



where the unit of electric field Ec{0) is defined in Eq. (4.20), and this will be a nontrivial test of the validity of 



the selfconsistent approximation applied to the TDGL equation and, following from it, the Boltzmann equation for 
fluctuation Cooper pairs. Now we are addressing cases important for the applications, in which the fluctuation 
superconductivity can be easily investigated. 

VI. STRIPED SUPERCONDUCTORS AND THICK FILMS 

The latest achievements in nanotechnology provide us in principle with a tool for a practical realization of so called 
striped superconductors with controlled parameters. For this it is necessary to cut stripes from a superconducting 
film using some appropriate lithographic technology. The amazing observation about the striped materials is that 
they are of "intermediate" dimensionality, i.e., they are neither 1-dimensional systems nor 2-dimensional ones. The 
closer the striped superconductor is to its critical temperature Tc, the more "perfectly" two-dimensional material 
it becomes, because the stripes become increasingly coherent. Analogously to the Lawrence-Doniach model for the 
layered superconductor we can describe this situation in terms of "dimensional crossover" . In the present section we 
show that the fluctuation conductivity longitudinal to the stripes can be derived following the standard procedure 
for layered superconductors. We have also to mention that probably some underdoped cuprates are naturally striped 
and this phenomenon has been at the center of the attention attracted by HTSC during the last few years. 

It is obvious from Section for example, that the LD conductivity for a layered superconductor can be natu- 
rally derived from a two-dimensional AL conductivity just by integrating that over the momentum in the direction 
perpendicular to the plane 

-..fe/) - /li-o (^ + 2^,/) - />»o (. + ^(1 - ™.)./) (6.) 

Thus the fluctuation conductivity for a striped superconductor reads 

o-stripcd(e,/) = ^^cTiD (e+ ^(1 - cos9)J^ (6.2) 

where s now stands for the period of the stripes. 

For a thick film with thickness dfiim we have to sum over the discrete spectrum of the energy associated with the 
motion in z-direction: 

^^yp^) = — — ^ ( — — ) ^z, (6-3) 

since 



2m*' ao V dR\ 



Pz = ^n,, 71, =0,1,2,3,.... (6.4) 

"film 

In the expression for the fluctuation conductivity of a thick film we have now summation instead of integration 



1 oo 

crfiim(e,/) = X!^2D ii + £z{Pz)/ao,f) = ^ ctsd (e+ «c(0)/(ifiim)^ n^, /) , (6.5) 

f^film ^ "film ^ 

which for zero electric field can be readily performed using 



a? + 1)^71^ 2a6tanh(^) 2a2 ' 

n=0 \ o / 



(6.6) 
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Finally we obtain the expression for the 3D conductivity of the film, which interpolates between 2D and 3D behaviour 
analogously to the conductivity for a layered superconductor: 

- '^-(^^ \^ + 2a.)tanh(|^)j ^ 1^ + ^) ^ ' ^'-'^ 

o-iD(e) = , = = 1 (6-8) 

Note that the thick film becomes "rightly" fwo-dimensional in the vicinity of Tc, thus the dimensionality of it decreases, 
whereas for the layered superconductor the dimensionality goes up from Z? = 2 to Z) = 3 as we approach the critical 
tenrnerature. In such a way a thick film of a strongly anisotropic layered superconductor, ^c(O) <C s <C dfiim, can have 
twotZI dimensional crossovers. 



VII. DETERMINATION OF THE LIFETIME CONSTANT tq 



In the current section we show that conductivity measurements in a strong electric field can serve as a method for 
probing fundamental properties of superconductors such as the lifetime constant of metastable Cooper pairs tq and 
the coherence length ^(0). We demonstrate that our theoretical results can be effectively used for experimental data 
processing and determination of both tq and ^(0). 

As a rule in an experiment the temperature dependence of a resistivity pcxp(T) is examined. The experimentally 
measured conductivity is consequently doxp = 1/Pcxp(2^)- The in-plane current, if we take into account the first 
nonlinear correction, can be written in the general form 



(7.1) 



since the nonlinear correction to the fiuctuation conductivity is Act/ — ~AE^. 

Of primary interest for us are the superconducting films. Let us consider the two-dimensional BCS-like supercon- 
ductor. In this case the parameter r that determines the effective dimensionality of the superconductor is zero, and 
the periodicity of the Lawrence-Doniach model s dfam (^fiim is the thickness of the superconducting film). Thus 



.(2D) 



^filmja: ■ 



the two-dimensional current is jx 

^From the expression for the LD fluctuation conductivity Eq (4.23) it follows: 



V(2D) 
Jx 



47r^ 



(8e)' 



Et 



(7.2) 



In order to study the fluctuation effect on conductivity one should plot first of alLthe paraconductivity contribution 
to the resistance l/(l/pcxp(T') — ^/pn{T)) as a function of T (see for examplell3, where indium/oxide films were 
examined). For InOx films the value of Tioi = 1.16 can be used as a tool for the determination of the in-plane 
coherence length ^ab(O). 

In the general case of a layered LD superconductor the coefficient in the nonlinear correction to fluctuation con- 
ductivity according to Eq. (4.23) reads 



9^2 J 5_„31 

8 " ^ 16 J 



[e{e + r)] 



7/2 



(7.3) 



In this case we have three coherence lengths ^a(O), ^^(O)? Cc(0): and the current anisotropy Jmax/ Jmin- 

In this paper we develop a model-free method for the determination of the lifetime of Cooper pairs tq and the 
coherence length from the experimental results for the fluctuation conductivity. The final result is derived on the 
basis of the Eq. (4.15) after three-fold integration over some time interval (ei, £2)- We obtain the following expression 
for the lifetime 



To 



Ci~e-e^f [-AafieJ)]de 



1/2 



\eEx\UO) a{e,) - \aie2) + (ei - e,) a'ie^) + U^i ~ e^f ^"(£2) 



(7.4) 
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In ( |7.4D we have the notations 



d 



a'{e) = -a{e), a"(e) ^ —a{e), 



de 



de 



(7.5) 



and £2 — 0.2 determines the upper bound on t he te mperature below which the fluctuation conductivity can be reliably 
measured. In order to simplify the formula (7^) we can make the approximation €2 oo. Since the fluctuation 
phenomena are negligible already for T ~Tc~ 0.15Tc the following asymptotic conditions can be imposed in the limit 

£2 OO, 



<j{e2) ^ 0, £2^7' (£2) ^ 0, ela"{e2) ^ 0. 



(7.6) 



Finally we obtain the formulae for the lifetime of metastable Cooper pairs and the in-plane coherence length, which 
can be applied for experimental data processing: 



To 



Trel 



1/2 



1 



U(o) = 

cf. Ref.i Eqs. (201-203). Here 



3 l&ksT 
2 ^|eS,|ea(0) \a{e 



1/4 



1/2 



(e-e) [-Aa„(£,/i)]d£ 



- Acr/i(£,/l) = — ^cr(£), £ 

is a nonlinear correction to the magnetoconductivity and the magnetic field 



is oriented perpendicular to the a6-plane, where 



-8,2(0) = -T, 



is the slope of the upper critical field and 



dB,2{T) 



dT 



$0 



2<,(0) 



$0 



2Trh 



(7.7) 
(7.8) 
(7.9) 

(7.10) 
(7.11) 

(7.12) 
(7.13) 



is the magnetic fiux quantum. 



VIII. CONDUCTIVITY CORRECTION BY DETECTION OF 3RD HARMONICS 



As an alternative method for probing the fundamental constants of the BCS theory we suggest the systematic 
investigation of the third harmonic of the electric field generated by a harmonic current. Third-harmonic measurements 
are easier to perform than those of resistivity, and, moreover, the effect arising from fiuctuations is exceptionally 
pronounced. 

In general AC current response and investigation of higher harmonics is a standard method for investigation of 
nonlinear effects on superconductivity. A homogeneous electric field 

E,{t) = Eocos{ijt), (8.1) 

for example, creates a small nonlinear response for the first harmonic and a cubic field dependence of the 3rd harmonic 
of the current 
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= jif cos{ujt) + jsf cos(3wt), 



where for the amphtudes of the harmonics we have 



1 



J3f 



tAE^ 



■2) 



.3) 



If necessary, a smooth analytical normal part of the nonlinear coefficient A can be subtracted from the experimental 
data vAexp 



ANiT)+Aie), An{T)^A + BT + CT^, 



(8.4) 



in order to extract pure fluctuation behaviour from the nonhncar coefficient of the conductivity correction 

Aaf = -A{e)El (8.5) 



If we use the so defined conductivity correction Aaf the electric field E^ is actually cancelled in Eq. (7.7) and we 
have to use the coefficient -4(e) in the expression for tq, i.e. 



To 



2 \e\UO) U(e) 



1 



{i-eyA{i)de 



1/2 



(8.6) 



There is no doubt that the electric field is a useful tool for a theoretical analysis but for the experimental realization 
of the suggested method we have to apply a harmonic current and to measure the harmonics of the voltage 



I{t) = Iq cos ujt, U [t) = Ui f cos Ljt + Us f cos 3ujt + C/5/ cos buit + 



■7) 



For small current amplitudes used to avoid heating of the sample the voltage response is in first approximation linear 
and we have Ohm's law for the resistance of the superconductor strip with length L, width w and thickness dfiim, 



Uif=R{T)Io, R{T)^p[T)- 



P{T) 



1 



Eo 



(T)^ 

Then the absence of the 3rd harmonic of the current, j^j = 0, according to Eq. 



L 



p{T)Io 
wdmm ' 



j) with Exit) = U{t)/L gives 



p{T)L 



to lowest order in /q, and finally we obtain 



A{e) 



U3f 



= 4- 



(8.9) 



(8.10) 



In this way the nonlinear coefficient necessary for the determination of the lifetime constant in Eq. ( |8.6| ) can be 
expressed through the electronically measured current Iq, voltage amplitudes U^f and Uif, and the geometrical 
parameters of the strip L, w and dfum- So the suggested experiment can be performed in every laboratory involved in 
investigations of superconductivity. 

Let us describe qualitatively the temperature dependence of the intensity of the 3rd harmonic when the temperature 
is increased. In the superconducting state the voltage response is negligible and the 3f signal will appear abruptly 
when we reach the critical temperature. After a sharp maximum at the 3f signal will decrease with smaller slope 
and in the normal region the 3f signal will be small again and created only by the 2f oscillations of the temperature and 
the temperature dependence of the resistivity. In short we predict that U3 / will have a A-shaped asymmetric critical 
singularity. The location of this A-point provides a new method for the determination of the critical temperature 
of superconductors based on the properties of fluctuation phenomena. Our selfconsistent theoretical calculation is 
applicable above Tc where C/3/ is much smaller than the value at the A-point, but still clearly detectable experimentally. 
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IX. DISCUSSION AND CONCLUSIONS 



Let us start analyzing the results derived with the momentum distribution Eqs. (2.11) and ( ^.13 ). We have a 
characteristic velocity related to the equilibrium distribution 



and Eq. ( 2.13 ) can now be rewritten as 



Tlx 

n{va;,e±) = - — ——2 — • (9.2) 

[vx/vc(e)\ +e + ej_/ao 



According to the Drude consideration, a small electric field creates a drift velocity 



^^drift(e) = = — (9.3) 

TO* m* e 

and the dimensionless electric field is just the ratio of those two velocities, 

Je = TT— > Vc(e) =i Vei'Fcrmiexp — 1. (9.4) 

^c\^) \ PFcrmi Impairing / 

^From a microscopic point of view the characteristic thermal velocity Wc(e) is proportional to the Fermi velocity 
times the small parameter of the BCS theory, the famous exponent which contains the density of states at the 
Fermi level and the matrix element of the pairing interaction. This order of magnitude estimation is applicable to 
anisotropic gaps as well. In addition we have a critical slowing down multiplier ^/e. Those two factors significantly 
decrease the characteristic velocity and make possible the experimental observation of the electric-field correction to 
the fluctuation conductivity proportional to ff. For normal metals the ratio t'drift/'^Fermi 1 is extremely small and 
only AC oscillations of the temperature T{t), mainly 2f, can create harmonics in the voltage response. This effect 
should also be carefully taken into account for cuprate films for which the thermal resistance between the substrate 
and the film can be very high; this will be the subject of another work. 

Let us also consider in short the tjr-quasiparticle criterion to check whether or not Cooper pairs are quasiparticles 
in the usual sense of condensed matter physics, i.e., Tp(tp — ijl) /Ti ^ \. For p = taking the microscopic value for 
the lifetime we have 

/ Na(e) (BCS) flo TT flo 1 
n h lb kb^c 16 riT 

This strong inequality means that fluctuation Cooper pairs are not quasiparticles. The notion of Cooper pairs is 
only a language to describe the properties of slowly decaying diffusion modes of the superconducting order parameter 
above Tc. Analogously the "mean free path" /(e) = w(e)r(e) is also very short relative to the correlation radius: 

l{e)_ ^ vjeMe) ^ h r(e) ^ ooTq ^ J_ -, ,q R^ 

^(e) e(e) m*^(e) ^{e) h 8 nr ^ ' ^ ' ^ 

As an illustration let us take a set of parameters corresponding to a high- Tc cuprate: Tc = 90 K, fc^Tc = 7.76 meV, 
TO*h = llTOe, where rUe is the mass of a free electron, ^ah(O) = 11 A. Then oo = ?i^/2TO*j^^j(0) = 2.86 meV and 
riT = ksTc/ao = 2.71 > 1. The last inequality ensures that n(e) = ny/e ^ 1 and this condition of applicability of 
Rayleigh- Jeans statistics justifies the treatment of the Ginzburg-Landau ^'-function as a classical complex field. 

In view of the new effects which can be predicted using the derived velocity distribution we consider the possibility 
of supercooling of the normal phase, cf. Ref.113, to be very interesting. In this case the fluctuation conductivity could 
create a negative differential conductivity, which opens perspectives for many technical applications. In order to 
prevent the nucleation of superconductivity from regions where the current densities and electric fields are very small, 
depairing impurities should be introduced in the contact area of the microbridge. This could be realized, for example, 
by evaporation of Ni on the wide area of the Al microbridge or by Mn ions in cuprate films. In both cases the central 
narrow region of the microbridge should be protected. In a sample prepared under these conditions the criterion of 
negative differential conductance can be easily satisfied and so the predicted generation of oscillations would probably 
be the best example of significant fluctuation effects in superconductors. 
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Whether the threshold for the generation regime can be described within the nonhnear theory or whether we need 
to calculate the fluctuation density and the renormalized temperature Tr = Tc(er + 1) in a selfconsistent way, depends 
on the numerical value of the Ginzburg number. In any case the analyzed solution of the Boltzmann equation suffices 
to predict a cross-over from positive to negative differential conductivity, as the field decreases. As a precursor of 
oscillations, when the electric field is decreased, due to strong nonohmic behaviour and low dissipation, the sample 
will be an excellent frequency mixer. The incipient bulk conductivity should always be taken into account, because 
the supercooled normal state is metastable and applying a voltage to the superconducting state leads as a rule to a 
space- and time-inhomogeneous phase. 

Some words should be added concerning the history of the kinetic equation introduced in 1876 by Boltzmann. This 
was the first use of probability concepts in a dynamical theory but the real recognition of the Boltzmann equation 
was stimulated by the electronic industry in the second half of the 20th century. Indeed, the Boltzmann equation is 
an essential tool for understanding how electronic devices work. 

In sum, in our paper we have used the Boltzmann equation not only as a didactical instrument but also as a means 
of deriving new results. The Boltzmann equation for fluctuation Cooper pairs is a consequence of the time-dependent 
Ginzburg-Landau equation. This equation has been derived from microscopic theory and in this sense the Boltzmann 
equation is a tool for the application of the microscopic theory of superconductivity. That is why the Boltzmann 
equation can be used to predict the results of new experiments, to help in their interpretation, and even to correct 
some previously obtained results derived in the frame-work of microscopic theory. 
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APPENDIX 

In the Boltzmann equation Eq. ( |2.5| ) the decay rate of fluctuation Cooper pairs depends on the kinetic energy £(p). 
As an illustration here we will analyze the isotropic GL model with 

1 1 e(p) -I- aoe , , , de{p) 

^ '^(P) = TT-T' P = ™v, v=— — . (9.7) 



t(p) To ao 2m* dp 

The Boltzmann equation is a dynamic equation and we consider it natural to trace the time evolution of the velocity 
distribution starting from some arbitrary initial distribution 

n(p,t = 0) =no(p). (9.8) 

As a rule dynamic problems can be more elegantly analyzed using canonical variables. That is why we will analyze 
the momentum distribution N{F,t) defined with respect to the canonical momentum 

P = p + e*A. (9.9) 

For a constant and space-homogeneous electric field E = const in the ip = gauge we have 

A{t)^-tE, A(i = 0) = 0, P(p,i) =p-e*Et, (^^^ ^ 1^^^, (^^^ = -e*E. (9.10) 

Now we will solve the Boltzmann equation using the distribution N{P,t). The physical quantity is the same - the 
number of Cooper pairs living at some momentum point - but mathematically the functions are different and their 
correspondence is given by 

N{P,t) = n{p,t). (9.11) 

After this change of the variables the sum of partial derivatives on the left-hand-side of the Boltzmann equation 
reduces to a usual derivative 
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d_ 

dt 



e*E- 



9p 



,./T,. N N dN dN dV dN dV 
^(P(p,t),0 = - + -.- + -.-.e*E 

dN{V,t) dNp{t) 



dt 



dt 



(9.12) 



The physical interpretation is very simple, in an external electromagnetic field the canonical momentum is conserved 
and reduces simply to a label with which the distribution can be parametrized. In this derivation it is essential to 
note that when a partial derivative is taken with respect to one argument, the other argument of the function is kept 
constant. As a consequence, {dn/dt)^ ^ {dN/dt)p. Taking into account that the decay rate and the equihbrium 
distribution are functions of the kinetic momentum p = P — e*A, the Boltzmann equation takes the form 



dNp{t) _ Np{t) - n(P - e*A) 



dt 



t(P - e*A) 



(p + e*my 

2m* 



aoe 



Npjt) ^ nr 
aoTo To 



Introducing now dimcnsionlcss time and dimensionless canonical momentum 

u=-, q = ^(0)P/;i, q = k-fu, 

To 

we find that in the ID case the Boltzmann equation takes the convenient form 

dNp{n) 



du 

One can easily check that the function 



[q + Ju) +e Np{u) + nT. 



(9.13) 



(9.14) 



(9.15) 



(g + fu'Y + eu' } du' 



no(q)exp<j {q + fuf 



satisfies this ordinary linear differential equation with the initial condition 

Nq{u = 0) = no{q). 
Now introducing a new dimensionless variable 



u = u — u 



(9.16) 

(9.17) 
(9.18) 



we find for the distribution nk{u) = Nq{u) with respect to the kinetic momentum k = q + fi 



nk{u) = riT 



exp 



(fc2 + e)u + kfu^ - -fu^ \ du 



+ no{k - fu) exp l-{k'^ + ()u + kfu^ - -fii^ 



(9.19) 



Note that the argument of no is a conserved quantity k — fu, cf. also Eq. (9.16). 

The transition to the limit i = tqu — > oo is very simple. After several relaxation times we have the stationary 
distribution 



n(k) = nk{u — > oo) — ut 



expj- (fc2 + e)u + fc/u2- i/2u3|dM. 



(9.20) 



Rest orin g k k^, e ^ e + w, and n(k) n{kx] w, /, e) we re turn to the solution of the static Boltzmann equation 
Eq. (3.4). Consequently the dummy parameter t = tqu in Eq. ( ^.llD has the meaning of the time interval between the 
birth of fluctuation Coop er pairs and the moment when we measure the current. Hence the integrand in the formula 
for the current Eq. (2.15) has the meaning of the part of the current given by Cooper pairs born a time span t before 
the moment of the measurement. We consider it interesting that according to Eqs. ( ^.18 , 5.19| ) this "age distribution" 
has a sharp J-like maximum at to = tqMo for the supercooled normal phase in small electric fields. 
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Finally, it is amusing to note that although the fluctuation Cooper pairs are neither particles nor quasi-particles 
the conservation of canonical momentum implies Newton's equation of motion 

^ = e*E. (9.21) 
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